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We derive expressions for the coupling coefficients for electromagnetic four-wave mixing in the 
non-linear quantum vacuum. An experimental setup for detection of elastic photon-photon scat- 
tering is suggested, where three incoming laser pulses collide and generate a fourth wave with a new 
frequency and direction of propagation. An expression for the number of scattered photons is derived 
and, using beam parameters for the Astra Gemini system at the Rutherford Appleton Laboratory, it 
is found that the signal can reach detectable levels. Problems with shot-to-shot reproducibility are 
reviewed, and the magnitude of the noise arising from competing scattering processes is estimated. 
It is found that detection of elastic photon-photon scattering may for the first time be achieved. 
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I. INTRODUCTION 

The linearity of Maxwell's vacuum equations does not 
allow electromagnetic waves propagating in vacuum to 
interact, and therefore photon-photon scattering is for- 
bidden in classical physics. However, according to the 
theory of quantum electrodynamics (QED), the quantum 
vacuum possesses nonlinear properties. Photon-photon 
scattering may therefore occur, owing to the interaction 
with virtual electron-positron pairs. It is thought that 
QED non-linear vacuum effects may be of importance in 
the neighborhood of strongly magnetized astrophysical 
systems, such as pulsars [y, Q and magnetar environ- 
ments 0, ■ With the rapidly growing power of present 
day laser systems 0, we are at the brink of making the 
non-linear properties of the quantum vacuum directly ac- 
cessible for observations Effects such as electron- 
positron pair creation and clastic photon photon scat- 
tering may even playan important role in future laser- 
plasma experiments Q . One such example would be laser 
self focusing 0, lldj | , where laser pulse compression gives 
rise to field strengths close to the critical field strength, 
E crit « lO^Vcm" 1 . 

Direct observation of elastic photon-photon scattering 
among real photons would, because of its fundamental 
importance to QED, be of great scientific importance. 
Throughout the last decades, several suggestions on how 
to detect elastic photon-photon scattering have been 
made. For example, using harmonic generation in an 
inhomogeneous magnetic field [Tlj , using resonant inter- 
action between eigenmodes in microwave cavities |l2lll3} . 
using ultra-intense fields occurring in plasma channels 
|l4j . as well as many others, see e.g. Rcfs. miii [13 ■ 
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However, so far no suggestions have led to actual de- 
tection of photon-photon scattering among real photons. 
Related processes, but physically different from elastic 
photon-photon scattering, is photon splitting (see 
also 0| ) and Delbriick scattering [2(j , of which the lat- 
ter has been detected using high-Z atomic targets |2Cj . 

The present paper is an expansion of Ref. Q , where the 
possibility to detect photon-photon scattering in vacuum 
using four-wave mixing is investigated and a concrete ex- 
perimental configuration is suggested where three collid- 
ing laser pulses stimulate emission in a fourth direction 
with a new frequency. Both theory, argumentation and 
the results are in the present paper presented in more 
depth than before, but the paper also features new as- 
pects, e.g. the analysis of the expected shot-to-shot re- 
producibility in the suggested experiment, the analysis of 
alternative experimental configurations and an extended 
discussion regarding noise sources. 

Using four-wave mixing to stimulate photon-photon 
scattering has the advantage of not being limited by the 
low scattering cross section, er 77 w 0.7 x 10 -65 cm 2 |2lj . 
in the optical range. The idea of QED four-wave mi xing 
was first studied by p^l . Further theoretical [2^, HI] 
and experimental [26j studies have been performed since 
then. The main obstruction for successful detection of 
scattering events has been the lack of sufficiently power- 
ful laser technology. Since present petawatt laser systems 
only operate at a single frequency (with harmonics avail- 
able through the use of frequency doubling crystals), the 
configuration of beams and the geometry of the setup also 
becomes important. We show that any two-dimensional 
(2D) setup is unlikely to produce a signal distinguishable 
from noise, for presently available power levels. This is 
in contrast to the 3D setup presented in and in this 
paper. Here we have calculated the coupling coefficient 
for four-wave mixing as a function of incident angles and 
laser polarization. This is done both for the general 2D 
setup and a specific 3D setup. We have also determined 
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an easy-to-use expression for the number of scattered 
photons with the 3D configuration, given specific beam 
data. The beam data parameters are chosen to fit the 
high-repetition rate Astra Gemini system (operational 
2007) located at the Central Laser Facility (Rutherford 
Appleton Laboratory), giving an estimated number of 
0.07 scattered photons per shot. Furthermore, a unique 
fingerprint of the QED process can be obtained through 
the polarization dependence of the number of scattered 
photons. 

In order to get a statistically sufficient number of scat- 
tered photons, a good shot-to-shot reproducibility is re- 
quired. The shot-to-shot reproducibility of the Astra 
Gemini system is analysed, and it is found that energy 
fluctuations will be small on a shot-to-shot basis. While 
vibrations could pose a problem, they are expected to be 
small due to the inherent structure of the Astra Gem- 
ini system and the experimental setup. Further ways 
of improving the vibration stability, if necessary, is also 
suggested. 

Due to a non-perfect vacuum in the interaction cham- 
ber, competing scattering processes such as Compton 
scattering and in principle collective plasma four-wave 
mixing will be present. An analysis, with the help of com- 
puter simulations, has been carried out in order to give 
a rough estimate of the signal to noise ratio. It is found 
that plasma cavitation caused by the strong laser pulses 
suppress the competing effects rather effectively close to 
the interaction region (see e.g. Ref. [13] )• Specifically, in 
our suggested experimental setup, this reduces the noise 
to a level approximately three orders of magnitude lower 
than the QED signal. Thus, we conclude that detection 
of photon-photon scattering will in principle be possible 
with the Astra Gemini system. 



On identifying an effective polarization and magneti- 
zation of the vacuum as 

P = -p [2(£ 2 -B 2 )E + 7(E-B)B] (2) 

and 



M= — \-2(E 2 -B 2 ) B + 7(E-B)El , (3) 

the equations of motion following from the Heisenberg- 
Euler Lagrangian can be written 

V • E = -4ttV • P, (4) 
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c dt 



V-B = 0, (6) 



1 <9E 1 dP 

V x B — = 4tt — — + 4?rV x M, (7) 

c dt c dt 

that is Maxwell's equations in the presence of a medium 
with polarization P and magnetization M. From those 
we can derive the effective wave equation 



□E = 4ttc 2 



V(V-P) 



Id fldP 



dt 



dt 



V x M 



> (8) 



where □ = df — c 2 V 2 . 

Next we consider three plane waves, representing the 
incoming laser pulses, with amplitudes allowed to have a 
weak space time dependence due to interactions, 



II. THEORY 

It is possible to effectively relate the properties of the 
QED-vacuum to those of a medium in ordinary classical 
electrodynamics. By integrating out all high energy de- 
grees of freedom, an effective field theory containing only 
the electromagnetic fields can be formulated in terms 
of the Heisenberg-Eulcr Lagrangian psl . l29l |. valid for 
field strengths below the QED critical field 10 16 Vcm -1 
and for wavelengths larger than the Compton wavelength 
10~ 10 cm J^J. The Heisenberg-Euler Lagrangian is given 
by 
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where £ = /ie 4 /457rm 4 c 7 , h = h/2n, h is Planck's con- 
stant, e is the magnitude of the electron charge, m is the 
electron mass, and c is the speed of light in vacuum. £ 
is the Lagrangian for classical electrodynamics, while SC 
represents the lower order non-linear QED-correction. 



Ej(r,t) = \ [E j (r,t)e iv ^ T - t) + E*(r,t) e"^^) , (9) 

where tpj = kj • r — ujjt and j = 1, 2, 3. Due to the cubic 
non-linearity of Q and J3J), we expect generation of a 
fourth wave with (104, h^) if the incoming wave vectors 
satisfy the matching conditions 



and 



ki + k 2 = k 3 + k 4 



UJ\ + uj 2 — UJ3 + UJ4, 



(10) 



(11) 



for some k 4 and 004 = clt4. In practice the interaction 
pulses will have a finite spectral width, Akj and equa- 
tions l |l and refers to the central peaks of the spec- 
tra. However, we assume that the spectral width is small 
(| Akj I <C |kj |), such that it can be incorporated into 
slowly varying amplitudes, \dEj/dt\ <C ujj \Ej \ , \VEj\ <C 
\kj \ \Ej\ and the derivatives can be taken to act only on 
the harmonic parts. The effect of a finite spectral width 
is to give a finite interaction time. Note that for a suf- 
ficiently broad spectra, more general methods must be 
used, see e.g. !3jj . From l|10|) and (|llf) we find that 
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only the resonant terms including the factors E1E2E3 
and E^E^Es are of importance, since all others will av- 
erage to zero over a short spacetime interval due to rapid 
oscillations. 

Hence, on neglecting the non-resonant terms, the wave 
equation JHJ) for the generated field E g takes the form 



□E 9 (r,t) 



^ulGE^Ele^-*- 



where 



E 3 (r,t) = -(E 3 (r,t) + E;(r,i)) 



(12) 



(13) 



and G is a geometric factor which only depends on the 
directions of the wave vectors and polarization vectors, 
and can be solved for from (JSJl . The driving of the initial 
waves can be neglected due to the weakness of the gener- 
ated field, and hence the energy in each pulse is constant 
during the interaction. 

The radiation zone solution to the wave equation i|12|) 
takes the form of a spherical outgoing wave multiplied 
by a direction dependent factor peaked in the resonant 
k4-dircction, 



E fl (r,i) 



Q e i(fc 4 r— Wit) 



< J (EiE 2 E 3 



e iki (^-*)' r 'dV', (14) 



where V' is the interaction region, and the amplitudes are 
to be evaluated at the retarded time 4r = t— |r — r'| /c. 

Equation l|14|) holds for any set of wave vectors sat- 
isfying (|10fl and l|ll[l . but for real experiments only a 
few configurations are of interest. In practice, the high- 
est power laser systems operates at a rather well-defined 
frequency, and by using frequency-doubling crystals we 
also have access to second harmonics, although with some 
power loss. As will be argued below, these restrictions in 
incoming frequencies make a 2D setup less competitive 
than a 3D setup. The 3D configuration achievable from 
a single laser through frequency doubling which is best 
suited for an experiment (allowed by the matching con- 
ditions i|lU|) and dJ) is given by the wavevectors 



ki = fcx 
k 2 = fcy 
ka = ¥ 

k4 = fcx + fcy — hki 



(15) 



and is illustrated in Fig. 2] 

The main advantages of using a configuration defined 
by l|15(l are that the wavelength and the direction of the 
generated wave are well separated from that of the in- 
coming laser beams. This is highly desirable since opti- 
cal components will always scatter and reflect some light, 
making it difficult to distinguish QED scattered photons 
of wavelengths close to that of the incoming laser beams. 
Furthermore, the geometry of the setup is not hindered 




FIG. 1: Configuration of the incoming laser beams and the 
direction of the scattered wave for the wave vectors defined in 
115H . which satisfies the matching conditions 1101 and 1111 . 



by the fact that the focused laser beams in reality are 
cone-shaped. Also, the polarizations can be chosen ar- 
bitrary since it is not necessary to prevent any counter 
propagating laser beam from entering back into the laser. 

In order to calculate the geometric factor G we define 
7j as the angle between the z-axis and the polarization 
vector of the Sj— field, in such way that it is positive in 
clockwise direction when looking in the kj -direction. (3$ 
is defined in the same manner but with respect to the 
x-axis. In the case of configurations confined to a plane, 
the wave vectors are given by 



kj = kj cos ( 



>jX + kj sine 



(16) 



where <pj is the angle from the x-axis to kj . Using these 
definitions, the general 2D and the s pec ific 3D geometric 
factors are given by (|A2|) and l|A3(l [32j . respectively, in 
the Appendix. \K2\ is valid for any set of wave vectors 
confined to a plane as long as (|10|) and (|11|) are satisfied, 
while i|A3() is valid for the specific set of wave vectors 
given by ||TBJ|. Note that the 2D case describes the general 
four-wave process in the center-of-mass frame of two 
incoming photons, with wave vectors and k 2 = — 
Thus, the general 3D geometric factor G can be found 
from the 2D one by applying Lorentz transformations. 

The number of scattered photons can be estimated 
from modeling the laser pulses as rectangular prisms with 
length L and quadratic cross section b 2 , inside which the 
field amplitudes are constant. Although not a completely 
accurate model of the pulse shape, it has the advantage 
of making the calculations transparent, and it will give a 
reasonable estimate of the number of generated photons. 

In the 3D case the interaction region will take the shape 
of a cube with side b, existing during a time L/c. By car- 
rying out the integration in (|14|l . the generated electric 
field can now be solved for, from which the generated 
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intensity and total power can be calculated. After mul- 
tiplying the generated power by the interaction time and 
dividing by the photon energy hio 4 , the estimated number 
of scattered photons per shot is found to be 
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(17) 



where G^d is the geometric factor |G3<j| which is given 
by (|A3|) . Pj is the power of the incoming pulses, A4 the 
generated wavelength, and with rf given by 
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(18) 



where I g is the generated intensity, I gtres its maximum 
value in the resonant direction, f{0,4>) = 1 — cos sin 
and g{0,4>) = sin sin 9. The approximate equality of 
if is accurate within 7% for 6/A4 < 30, i.e. 6/A3 < 10 
where A3 is the fundamental wavelength. 77 corresponds 
to roughly half the angular width of the central interfer- 
ence peak of the generated intensity. 

The scaling of the number of generated photons with 
respect to G 3 d is decoupled from the pulse model. Hence 
it is possible to test the theoretical predictions by vary- 
ing the polarizations of the incoming pulses. However, 
changing the linear polarization of one of the laser beams 
could create a false periodic signal, as the polarization af- 
fects the propagation direction of the competing Comp- 
ton scattered photons. A low noise contribution is there- 
fore a necessity. The upper panel of Fig.|2shows how the 
scattering number depends on the polarization angle of 
wave three, with the polarization of the other two waves 
fixed at 71 = and 72 = 7r/2. An optimal choice of po- 
larization angles is given by 71 = and 72 = /?3 — 7r/2, 
for which G\ d = 0.77. 

Should the polarization dependence from experiments 
deviate from that predicted by theory, it would be in- 
teresting to know whether this deviation also suggests 
deviation from Lorentz invariance. Lagrangians of Born- 
Infeld type are Lorentz invariant and can give us some 
insight in how to interpret a deviation. Nonlinear elec- 
trodynamics was suggested by Born & Infeld [23, and 
later reformulated by Dirac j3J], already in 1934, as a 
means for obtaining finite field energies around singular 
charges. In essence, the Born-Infcld Lagrangian builds 
on the Lorentz invariants constructed from the Maxwell 
field tensor, and can be written [33 
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where k is the relevant coupling constant for the modified 
electrodynamics. The Born-Infeld form of Lagrangian 
(|19J) also occurs as the effective field limit of quantized 
strings [3(| (see Ref. [37J for a review). In such effec- 
tive Lagrangians the field strength coefficients contain 
the string tension. We note that l|19|) gives no birefrin- 
gence in static fields, a unique property of the Born- 
Infeld Lagrangian j35j. In general one speaks of Born- 
Infeld type Lagrangians when we have Lagrangian built 
from the Lorentz field invariants, and the coefficients are 
kept arbitrary, see (|A1|) . A Lagrangian of this form af- 
fects the Geometric factor, and thus the polarization de- 
pendence of the signal as described in Appendix A. Any 
detected signal that can be fit to follow the polarization 
dependence of the geometric factors for some X and Y 
is thus consistent with Lorentz invariance. 



III. EXPERIMENTAL PARAMETERS 

The Astra Gemini Laser will generate two indepen- 
dently configurable 0.5 PW laser beams of wavelength 
800 nm. Each pulse will contain a total energy of 15 J, 
and focused intensities over 10 22 Wcm -2 will be reached. 
The repetition rate is expected to be three shots per 
minute. Using these values, the spatial dimensions of 
the pulse model are taken as b = 1.6 /im and L = 10 fim., 
which gives rf = 0.025. The 3D configuration described 
by (|15fl can be achieved if one of the laser beams is fre- 
quency doubled and split into two beams. The estimated 
loss of power when frequency doubling is approximately 
60%, and hence the power of the incoming beams are 
given by P l = P 2 = 0.1 PW, P 3 = 0.5 PW. The number 
of QED scattered photons, using optimal polarizations, 
would then be N$d = 0.07 and their wavelengths will be 
centered around 267 nm. The lower panel of Fig. [21 and 
Table [I] shows how the scattering number increases with 
increasing laser power, when the focused beam width is 
kept constant. 

We note that in order to get a statistically sufficient 
number of scattered photons, the conditions must be re- 
produced from shot to shot. The Astra Gemini system is 
designed to operate in saturation, so energy fluctuations 
are likely to be less than 5% on a shot to shot basis. 
As for the pulse length, it primarily derives from static 
and fixed effects and fluctuations are expected to be less 
than 5%, giving a maximum total fluctuation in power 
of about 10%. The most important effect is that of spot 
and pointing stability. Through the use of adaptive op- 
tics, there will be active control of the spot size, shape 
and to some limited extent position. Controlling the rel- 
ative pointing of the three beams after they have been 
split and frequency converted is a challenge. This will be 
done at large aperture so there is an inherent sensitivity 
to small pointing variations arising from vibration. Vi- 
bration could pose a problem, but the experiment should 
be designed in such a way that any given mode of vi- 
bration produces a similar far field displacement, i.e. all 
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TABLE I: The table shows more precise values for the num- 
ber of scattered photons in Fig. 2, for some different laser 
powers. N a is the number of scattered photons for a system 
where two source beams are used and Ni is the number when 
three source beams are used. 0.066 is the expected number of 
scattered photons per shot with the Astra Gemini system. 
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2.7 x 10 5 



beams move in the same direction vertically or horizon- 
tally against a stimulus. Furthermore, the Astra Gemini 
system has been designed on very thick concrete slabs so 
vibration after this split is likely to be a minimum prob- 
lem. Should this prove not to be the case, then com- 
mercially available stabilising systems would be used to 
stabilise the pointing of the beams. Finally, all of these 
parameters will be measured on each shot in order to sort 
the data accordingly. Below we will show that already 
the scattering number N^d = 0.07 for the AG system will 
exceed the estimated noise level. 




Source beam power [PW] 



FIG. 2: The upper panel shows the number of scattered pho- 
tons N3d per shot, normalized by the number of photons iV ma x 
per shot for an optimal choice of polarization, as a function of 
polarization angle fig of the wave with direction kV The lower 
panel shows -/V max predicted by l|17jl when increasing the laser 
power while keeping the beam width constant at 6 = 1.6 /im, 
for a system where two source beams are used and one of 
the beams is split into two (solid line) and when three source 
beams are used, hence no beam splitting is required (dashed 
line). 



IV. ALTERNATIVE SETUPS 

For comparison, let us consider a 2D configuration 
where two parallel beams collide head on with a third 
beam, all of them with the same wavelength. The inter- 



action region will then take the shape of a rectangular 
prism of width b and, with respect to the retarded time, 
length L/2, and exist for some time interval L/c. Since 
the pulse length often greatly exceeds the pulse width, 
this configuration gives an optimal interaction region and 
also generates the most scattered photons, TV = 2.4, for 
an optimal choice of polarization angles (71 = 73 = and 
72 = tt/2). However, the scattered photons are emitted 
along the beam axis and have the same wavelength as 
that of the incoming beams, making detection impossible 
in practice. The directions of the incoming and scattered 
waves can in principle be separated if their frequencies 
are shifted. The large cone angle of the laser beams and 
the need of having a well separated scattered frequency 
however requires a large deviation from the original con- 
figuration, thus destroying its nice features and conse- 
quently bringing down the scattering number. This also 
requires a greater change in frequency of the incoming 
beams, which may cause a problem with finding suitable 
lasers. 

There is, however, a possible 2D configuration where 
only the fundamental frequency and its second harmonic 
is used, with the wavevectors 
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k 2 = 
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(20) 



as illustrated in Fig. EI The geometric factor for this con- 
figuration becomes approximately G\ d = 5.5, for an opti- 
mal choice of polarizations angles (7! = 72 = 73 = 3n/A). 
Although this looks superior to the 3D setup, there are 
some practical problems. To extract the QED photons 
in this 2D setup, a beam splitter needs to be introduced 
in the path of k 3 , before the interaction region, bringing 
down the scattering number with a factor 0.25. Also, the 
interaction region in this case becomes somewhat smaller 
than a cubic one, further reducing the scattering number. 
Consequently, the estimated detectable signal will be of 
the same order as that of the 3D setup. The polarization 
dependence of the signal with this configuration is, how- 
ever, much more pronounced than for the 3D one and 
it nearly vanishes at minimum. This setup could be an 
alternative to the 3D configuration presented above, but 
a major remaining problem is to ensure that none of the 
higher harmonic photons of the fundamental laser beam 
is reflected by optics back to the detector, drowning the 
signal in noise. 

There are other possible experimental configurations, 
both 2D and 3D, that we have not investigated in this pa- 
per. It is, for instance, possible to use third harmonics or 
higher, to construct setups with the scattering direction 
distinct from the incoming photon directions. However, 
this requires additional frequency doubling crystals re- 
ducing the power, and we expect these configurations to 
be less competitive. It is also possible to use configura- 
tions consisting of two or more lasers of different frequcn- 
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FIG. 3: An alternative experimental configuration for detec- 
tion of photon-photon scattering through four-wave coupling. 
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cies to generate QED scattered photons with wavelengths 
well separated from the harmonics of the incident beams, 
as used in the 3D configuration in Ref . [2i| . Since present 
laser systems have sufficient power only at single frequen- 
cies (with harmonics available using frequency-doubling 
crystals) we expect these configurations to be less com- 
petitive. Hence, the 3D configuration, where just the 
second harmonic is needed and the scattering direction 
is distinct from the directions of the incoming photons, 
seems to be the optimal setup from an experimental point 
of view. 



FIG. 4: The maximum intensity felt by the electron relative 
the peak intensity of the laser pulse, for different laser inten- 
sities. No initial radial displacement of the electron from the 
propagation axis of the laser pulse. Laser wavelengths are 
800 nm and 400 nm. 
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Since a perfect vacuum in the interaction chamber is 
in practice impossible to achieve, there will be competing 
scattering processes present. The effect of these processes 
will rather effectively be suppressed due to plasma cavi- 
tation caused by the strong laser pulses. Below, with the 
help of computer simulations, we give a rough estimate 
of the magnitude of the competing scattering events re- 
sulting from Compton scattering and collective plasma 
four-wave mixing. This is in order to determine whether 
it is theoretically possible to distinguish the QED scat- 
tered photons from noise. 
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FIG. 5: The maximum intensity felt by the electron relative 
the peak intensity of the laser pulse, for different initial radial 
displacements from the axis of propagation of the laser pulse. 
The intensity is 5 x 10 21 Wcm -2 , and the wavelengths are 
800 nm and 400 nm. 



A. Ponderomotive force 

To investigate the effect of the ponderomotive force 
near the focal spot, we have given the laser pulse a Gaus- 
sian shape. With the beam parameters given for the 
Astra Gemini system, simulations show that an electron 
exposed to such a laser pulse near the interaction region 
will be forced out of the central beam already during the 
initial part of the pulse. Fig. 0] and [5] show the maximum 
intensity felt by the electron relative to the peak inten- 
sity of the laser pulse, for laser pulses with wavelengths 
800 nm and 400 nm. In Fig. the intensity of the laser 
pulse is varied and the electron has no initial radial dis- 
placement from the propagation axis of the laser pulse. 
In Fig. [S] the initial radial displacement is varied for an 
intensity of 5 x 10 21 Wcm -2 . 



Even though the electron blow out leaves a much 
denser ion plasma, the relative contribution to Comp- 
ton scattering from ions compared to electrons scales as 
the square of the mass ratio. We therefore neglect the 
ion contribution to the competing noise. 



B. Compton scattering 

For the intensities considered in this experimental pro- 
posal, both linear and non-linear Compton scattering 
will be present. The non-linear Compton scattering 
for high intensities have been described both classically 
|38ll39j | and semiclassicaly |40j for an infinite plane wave. 
The theory predicts generation of a scattered spectra 
with sharp peaks at each harmonics. This could pose 
a major problem since 3rd order Compton scattering 
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of the fundamental frequency generates photons with a 
wavelength that can be confused with the QED signal. 
However, with a focused laser pulse of this intensity, the 
broadening of the spectra will be tremendous. This is 
confirmed in experiment by |4lj . where the peaks at the 
harmonics are not distinct. 

We are only interested in the number of Compton scat- 
tered photons with energies in the same range as the 
QED scattered photons. Non-linear Compton scatter- 
ing of 3rd order and higher for 400 nm beam and 4th 
order and higher for the 800 nm beam are unlikely to 
yield photons with energies possible to confuse with the 
QED ones, and can therefore be disregarded. Non-linear 
Compton effects does not become significant until the 
dimcnsionlcss parameter, 



scattered photons, e sc , will then be within the range 



eE r 



ujmc 



(21) 



approaches or exceeds unity (21J, |42j . Here E rms is the 
root-mean-square value of the electric field amplitude of 
the laser pulse. Near the interaction region, electrons 
will be forced out of the central beam already by the 
leading edge of the pulse, thus experiencing an intensity 
significantly less than the maximum intensity of the laser 
pulse. With the given beam parameters and a maximum 
intensity I max of 5 x 10 21 Wcm -2 , most electrons will ex- 
perience an intensity less than 0.5% of J max (see Fig. [SJ, 
giving rj w 1. For most electrons close to the interac- 
tion region, 2nd and 3rd order processes will give a total 
cross section of the same order of magnitude as that of 
1st order Compton scattering [3^.l43l| . These higher order 
Compton scattering processes may be the largest source 
of noise at wavelenths close that of the QED photons. 
However, we note that we only need a order of magni- 
tude estimate of the noise level. Because of that, and 
because of the complexity of the problem we are satisfied 
with estimating the contribution from 1st order compton 
scattering. 

From now on we will refer to 1st order Compton scat- 
tering simply as Compton scattering. We use a hybrid 
model, where a classical electromagnetic wave accelerates 
an electron which in its restframe Thomson scatters pho- 
tons. This opens for inverse Compton scattering, where 
net energy is transferred from the electron to the pho- 
ton. We assume that we can filter out and detect only 
the photons within the wavelength interval ±50 nm cen- 
tered at the wavelength of the QED scattered photons, 
267 nm. We calculate the scattering cross section for in- 
verse Compton scattering into this wavelength interval, 
a a, in order to give a rough estimate of the noise caused 
by Compton scattered photons. For simplicity, the inci- 
dent light is assumed to be monoenergetic with energy e. 
We further assume that in the rest frame of the electron, 
the scattering is elastic and isotropic Thomson scatter- 
ing, which is a decent approximation as long as je <C mc 2 
|44j . where 7 is the relativistic factor. The energy of the 



1-/3 esc 1 + 
1 + 13 e 1-/3' 



(22) 



where (3 = u/c and u is the electron velocity. The cross 
section, <7a, for inverse Compton scattering into the en- 
ergy interval defined by e iow and e up , e < e iow < e up , can 
be shown to be (see [44)1 

(JA = (23) 

if (l+/3)/(l-/3) < e low Je, 

o T (e (1 + (i? h ,„ v , n m 4™ 



if eiotv/e < (1 + /?)/(! - /3) < e up /e, and 



(24) 



<7T 



4e 7 2 /3 2 



(1 + 



f 2 -e 2 



low 



2c 



(25) 

if e up /e < (1 + f3) /{I - [3). Here a T = 6.65 x 10" 25 cm 2 
is the Thomson scattering cross section. 

Using the equations of motion for a free electron in an 
electromagnetic field, dp a /dr = — (e/m) x F a bp b , the mo- 
tion of the electron can be simulated and by using l|23|) 
to (I25fl the average cross section during a wave period 
for scattering into a given wavelength interval for a given 
field intensity can be estimated. Here p a is the electron 
four-momentum, r is the electron eigentime, and F a b is 
the electromagnetic field tensor. The average scattering 
cross section for scattering of incident light with wave- 
length 800 nm and 400 nm into the wavelength interval 
216 — 316 nm for different field intensities is plotted in 
Fig. El 

An alternative model to the Compton model has also 
been used to verify the magnitude of the competing scat- 
tering events. The model is based on the principle that 
the laser pulse gives rise to electron oscillations, and for 
highly relativistic electrons, the parallel component of the 
acceleration can be neglected compared to the orthogo- 
nal one. The radiation is then synchrotron radiation for 
a particle instantaneously moving in a circular motion of 
radii p = u 2 /iij_ [45| . Simulations have shown that the 
scattering cross section with this model is roughly the 
same as that from the Compton model. 



C. Collective plasma effects 

In addition to the QED non-linearities, medium non- 
linearities can also contribute to scattering of photons. 
For example, in a similar configuration as that consid- 
ered by us, Ref. j2(| has detected photons scattered due 
to third order susceptibilities in a neutral nitrogen gas. 



8 



x1CT ! " 




.— 


■ X = 800 nm 


6 




• 


. X = 400 nm 






• 




5 




• 

• 








• 




4 












• ■ • 




3 




■ 




2 




■ ■ 








• 




1 






• 

• 








■ . 


16 


17 18 19 


20 






log of Laser Intensity [Wcm 


- 2 ] 



FIG. 6: The average cross section, <ta, for inverse Compton 
scattering of incident light with wavelength 800 nm and 400 
nm into the wavelength interval 216 — 316 nm, as a function 
of increasing intensity of the incident pulse. 



In our case, with field strengths well above the ioniza- 
tion threshold, the pulses will interact with a plasma 
rather than a gas, and in principle plasma collective non- 
linearities could compete with the QED effect. However, 
with a sufficiently low pressure (of the order 10 _9 torr), 
combined with the hole formation caused by the pondcro- 
motivc force, only a few electrons will be present within 
the pulsc(s), and none within the central parts of the 
focused pulse. As a consequence, collective plasma phe- 
nomena will not contribute significantly to the scattering 
made by the electrons. 



D. Estimate of Compton scattered photons 
scattered into the detector 

We now have the tools to make a rough estimate of 
the number of Compton scattered photons. If the elec- 
trons were distributed homogeneously within a pulse of 
N p photons, the number of scattered photons during a 
time r would be N sc = an e cN p T. However, because of 
the ponderomotive force, most of the photons will not 
be exposed to any electrons in the region near the focal 
point, where the intensity is the greatest. We thus in- 
troduce a parameter k, which is the fraction of photons 
that do experience electrons. With a maximum intensity 
of 10 22 Wcm -2 in the interaction region, we find from 
Fig. 0] and [S] that the effect of the ponderomotive force 
is crucial. We estimate that more than 99% of the beam 
photons do not encounter any electrons, and we reduce 
the scattering rate correspondingly, i.e. k = 0.01. 

Only a fraction of the scattered photons will have en- 
ergies possible to confuse with QED scattered photons, 
and within that energy range we introduce the subscript 
A on the cross section a. Noting that due to plasma cav- 
itation, the intensity felt by the electrons is significantly 
less than the peak intensity, we find from Fig.[S]the aver- 
age cross section <ta ~ 10 -26 cm 2 . From Fig. 0] and [5] we 
find that the effect of the ponderomotive force is much 
more crusial for the 800 nm pulse than for the 400 nm 
one. Furthermore, the average cross section a a is smaller 



for the 800 nm photons (sec Fig. EI) • Guided by this, we 
neglect the scattering from the 800 nm beam compared 
to that from the two 400 nm beams. The number of 
photons in one beam then becomes about N p w 6 x 10 18 
(40% of 7.5 J at 400 nm). 

Of the scattered photons, only a fraction, 6, will be 
scattered into the detector. Electrons accelerated by 
the laser pulse will oscillate relativisticly in the plane 
spanned by the polarization vector and the wave vector. 
Because of the relativistic beaming effect, the Compton 
scattered photons will be confined close to this plane. If 
we assume that the detector sees all photons from the 
QED scattered wave within a cone angle of 30 degrees, 
this would yield 5 = 0.017 if all photons were scattered 
isotropically. However, when considering that electrons 
are highly relativistic near the interaction region and tak- 
ing into account the relativistic beaming effect, we esti- 
mate S ~ 0.01. The expression for the number of Comp- 
ton scattered photons per beam reaching the detector is 
then 

Ndet = n5<j/±n e cN p T. (26) 

We assume that the interaction region is imaged onto 
the detector in such a way that the detector only sees 
a part of the interaction region about 10 /jm in di- 
ameter. This gives us an interaction time of about 
t ~ 3 x 10 -14 s, during which competing Compton pho- 
tons can be scattered into the detector. We further as- 
sume an unperturbed vacuum around 10~ 9 torr, giving 
n e m 5 x 10 9 cm" 3 . The number of Compton scattered 
photons that the detector might confuse with QED scat- 
tered photons, from both of the 400 nm beams, is then 
about Ndet ~ 5 x 10 ~ 5 , about three orders of magni- 
tude lower than the expected number of QED scattered 
photons. Here the contribution of higher order Compton 
scattering has not been included. Even if the non-linear 
contribution to the noise exceeds the linear one with an 
order of magnitude, it is still well below the QED signal. 
Should the noise level in a real experiment still prove to 
be too severe, it can be further reduced by ionizing the 
interaction region with a strong laser pulse and extract 
the electrons with a static electric field shortly before 
performing the QED experiment. 

VI. CONCLUSIONS 

In this paper, we have investigated the possibility to 
detect QED elastic photon-photon scattering through 
four-wave mixing. In the suggested 3D experimental 
setup, three colliding laser pulses stimulate emission of 
a fourth electromagnetic wave with a new frequency and 
a new direction of propagation. Considering that present 
day petawatt laser systems only operate at a single fre- 
quency (with harmonics available through frequency dou- 
bling crystals), it has been shown that a 3D configura- 
tion will be able to produce a measurable signal for the 
Astra Gemini system and similar high power systems. 
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Moreover, by changing the polarization of one of the in- 
coming beams, a unique fingerprint characterizing elas- 
tic photon-photon scattering can be obtained. For given 
data of the laser beams, together with the coupling co- 
efficients, an estimated number of scattered photons was 
obtained. With the parameters chosen to fit the Astra 
Gemini system at the Central Laser Facility (Rutherford 
Appleton Laboratory), the estimated number of scat- 
tered photons per shot was 0.07. Two main obstacles 
for successful detection were identified, the shot-to-shot 
reproducibility and the noise from competing scattering 
events. We have argued that pulse energy fluctuations 
will be low, and that vibrations are likely to be a small 
problem due to the inherent structure of the Astra Gem- 
ini system. Means of further improvements to reduce vi- 
brations have also been suggested. Due to a non-perfect 
vacuum in the interaction chamber, there will be compet- 
ing scattering processes such as Compton scattering and 
collective plasma four-wave mixing present. Through 
the use of computer simulations, it was found that the 
competing effects arc suppressed rather effectively due 
to plasma cavitation caused by the strong laser pulses. 
With the use of ultra-high vacuum technology, it was 
found that the noise from competing scattering processes 
would roughly be about three orders of magnitude lower 
than the QED signal. Means of further reducing the 
noise has also been suggested. We conclude that detec- 
tion of elastic photon-photon scattering will in principle 
be possible with the Astra Gemini system operational in 
2007. 
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APPENDIX A: THE GEOMETRIC FACTOR 

The explicit expressions for the general 2D and and the 
specific 3D geometric factors, G2d and G$d respectively, 
is found below for a Born-Infeld type Lagrangian 



1 

8^ 



(E 2 -B 2 ) 



8tt 



X (E 2 - B 2 ) 2 + Y (E • B) : 



_ (Al) 

where A and Y arc arbitrary coefficients. For the 
Hcisenbcrg-Euler Lagrangian, X — 1 and Y = 7. 



1 f \Y 
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